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Abstract 


Quantitative  Feedback  Theory  and  the  improved  method  Quantitative  Feedback 
Theory  are  enhanced  to  include  the  rejection  of  disturbance  at  the  system  output.  The 
enhanced  Quantitative  Feedback  Theory  and  improved  method  Quantitative  Feedback 
Theory  processes  are  applied  to  the  design  of  an  automatic  flight  control  system  to 
regulate  position  of  the  C-  135B  fuel  receiving  aircraft  relative  to  the  tanker  during  air-to- 
air  refueling.  A  simple  feedback  control  system  is  developed  that  will  achieve  stable 
position  regulation.  State-space  aircraft  models  are  generated.  An  "inner  loop"  autopilot 
is  designed  to  reduce  the  piant  cutoff  frequency  and  provide  the  system  inputs  for  the 
Quantitative  Feedback  Theory  compensators.  Disturbance  models  representing  disturbance 
due  to  wind  gusts  and  refueling  are  developed.  The  flight  control  system  is  designed 
using  the  enhanced  Quantitative  Feedback  Theory  equations.  Linear  simulations  are 
performed  on  MATRIXx,  and  nonlinear  simulations  are  run  on  EASY5x.  The  results  of 
the  simulations  show  excellent  results.  The  simulation  results  indicate  that  air-to-air 
automatic  flight  control  system  are  technically  achievable,  and  that  implementation  in 
US/  F  aircraft  is  possible. 
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DEVELOPMENT  OF  AN  AIR-TO-AIR  REFUELING 
AUTOMATIC  FLIGHT  CONTROL  SYSTEM  USING 
QUANTITATIVE  FEEDBACK  THEORY 

/.  Introduction 

1.1  Background 

The  United  States  Air  Force  (USAF)  has  the  mission  of  transporting  large  quantities 
of  material  and  troops  over  great  distances.  To  meet  this  requirement,  the  USAF 
maintains  a  fleet  of  large  cargo/transport  aircraft.  Refueling  these  aircraft  during  flight 
provides  unlimited  range  of  operation  for  this  fleet  of  aircraft.  However,  long  flights  and 
multiple  air-to-air  refuelings  can  seriously  strain  and  fatigue  the  pilot,  decreasing  flight 
safety,  and  extending  recovery  time  between  missions. 

Cargo/transport  aircraft  are  generally  large  and  have  high  moments  of  inertia. 
Piloting  a  large,  high  inertia  aircraft  during  air-to-air  refueling  requires  intense 
concentration.  The  pilot  must  maintain  a  very  precise  position  relative  to  the  tanker. 
He/she  maintains  position  visually,  applying  the  appropriate  control  inputs  when  changes 
in  position  occur.  The  pilot  must  compensate  for  changes  in  aircraft  dynamics  due  to 
talcing  on  fuel,  specifically,  changes  in  center-of-gravity  and  moments  of  inertia  I„  and 
I„.  Besides  dynamic  changes,  the  pilot  must  contend  with  maintaining  position  in  the 
presence  of  wind  gusts.  Since  these  aircraft  can  take  on  large  amounts  of  fuel,  up  to 
250,000  pounds,  air-to-air  refueling  can  take  up  to  30  minutes.  Compound  this  over  long 
flights  and  multiple  refuelings,  the  pilot's  fatigue  level  increases  and  can  reach  an  unsafe 
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level,  endangering  the  flight  crew,  and  possibly  impact  USAFs  capability  to  perform  its 
mission. 

One  way  to  ease  the  pilot  workload  is  to  implement  an  automatic  flight  control 
system  (AFCS)  for  air-to-air  refueling.  The  AFCS  needs  to  be  able  to  maintain  a  precise 
position  of  the  fuel  receiving  aircraft  (receiver)  relative  to  the  tanker  in  the  presence  of 
such  disturbances  as  wind  gusts  and  mass  changes.  In  this  thesis,  an  AFCS  is  designed 
to  precisely  regulate  position  relative  to  the  tanker  by  applying  the  control  synthesis 
techniques  of  Quantitative  Feedback  Theory  (QFT).  QFT  is  a  design  technique  developed 
by  Dr.  Isaac  Horowitz  that  quantifies  the  plant  uncertainties  and/or  disturbances  and  uses 
this  information  to  design  feedback  compensation  to  achieve  stability,  robustness,  and 
desired  system  performance  in  the  face  of  parametric  uncertainty  and  disturbances [2]. 

For  the  regulation  problem,  a  multiple-input  multiple-output  (MEMO)  QFT  design 
method  is  developed  to  address  the  disturbances  entering  die  system  at  the  output.  Dr. 
Horowitz's  original  MIMO  QFT  development  modeled  all  disturbances  as  entering  the 
system  at  the  input  of  the  plant  [3].  In  this  thesis,  MEMO  QFT  is  expanded  to  account 
for  disturbances  at  the  system  output. 

The  sponsor  for  this  thesis  is  the  Flight  Dynamics  Directorate,  Wright  Laboratory, 
Wright-Patterson  AFB,  OH. 

1 .2  Problem  Statement 

During  air-to-air  refueling,  the  receiver  aircraft  will  change  position  relative  to  the 
tanker.  The  pilot  must  pay  close  attention  and  take  corrective  action  to  maintain  position. 
Excessive  changes  in  position  will  disconnect  the  boom  from  the  receiver,  An  AFCS 
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must  be  designed  to  regulate  the  receiver's  position,  thus  reducing  the  pilot  workload  and 
fatigue  factor.  By  using  MIMO  QFT,  an  AFCS  is  designed  that  operates  throughout  the 
range  of  the  changing  aircraft  dynamics  and  rejects  disturbances  including  those  at  the. 
output  [3]. 

13  Assumptions 

The  following  assumptions  are  made  for  this  thesis. 

•  Only  the  desired  outputs  are  of  interest  for  final  performance. 

•  The  robustness  of  the  CAD  packages  used,  MIMO/QFT,  EASY5,  MATRIX,  and 
Mathematica  are  adequate  for  the  design  process. 

•  Position  of  the  receiver  aircraft  relative  to  the  tanker  during  air-to-air  refueling  can 
be  accurately  measured. 

The  first  assumption  is  required  in  applying  MIMO  QFT.  The  second  assumption 
is  concerned  with  the  limits  of  CAD  packages  and  their  numerical  robustness.  The  third 
assumption  is  required  because  no  sensors  are  currently  in  place  to  measure  the  position 
of  the  receiver  relative  to  the  tanker. 

1 .4  Research  Objectives 

This  thesis  (1)  develops  aircraft  models  for  the  QFT  design  from  a  published 
document  [1]  containing  C-135  cargo  aircraft  stability  derivatives  tables  and  plots,  (2) 
presents  a  design  for  multi-channel  control  laws  using  MEMO  QFT  for  several  flight 
conditions  with  special  emphasis  on  changing  aircraft  center-of-gravity  and  weight,  (3) 
simulates  the  design  for  linear  and  nonlinear  performance  on  MATRIX,,  and  nonlinear 
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performance  on  EASY5x,  (4)  evaluates  the  new  control  law,  and  (5)  validates  the 
development  of  MIMO  QFT  design  with  disturbances  at  the  output. 

7.5  Scope 

This  thesis  applies  the  MIMO  QFT  technique  to  the  design  of  an  AFCS  regulator  for 
the  three-dimensional  separation  (x,  y,  and  z)  of  a  C-135B  (receiver)  cargo  aircraft 
relative  to  a  tanker.  The  AFCS  applies  only  to  the  receiver  and  is  independent  of  the 
tanker  in  as  much  as  the  tanker  is  used  as  the  point  of  reference.  MIMO  QFT  design 
techniques  are  developed  for  disturbances  at  the  system  output.  The  present  MEMO/QFT 
CAD  package  [4]  is  modified  to  incorporate  the  output  disturbance  development.  The 
AFCS  is  designed  to  reject  disturbances  at  the  x,  y,  and  z  outputs  in  order  to  keep  the 
receiver  aircraft  in  a  volume  specified  as  the  area  of  boom  operation  [1],  Models  are 
developed  for  disturbances  due  to  wind  gusts  and  receiving  fuel.  The  MIMO  QFT  plant 
is  the  C-135B  bare-aircraft  model  augmented  by  a  typical  Mach-hold,  altitude-hold,  wing- 
leveler  autopilot.  QFT  compensators  control  the  reference  signal  of  the  autopilot  to 
achieve  "formation  maintenance"  during  air-to-air  refueling.  The  control  system  is 
simulated  for  linear  performance  in  MATRIX,.  A  full  six-degree-of-freedom  nonlinear 
simulation  is  performed  in  EASY5x. 

1.6  Methodology 

The  design  approach  requires  six  steps: 

•  Generate  linear  time-invariant  (LTI)  state-space  models  of  the  C-135  for  different 
weights  and  center  of  gravity. 
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•  Implement  a  Mach-hold,  altitude-hold,  wing  leveler  autopilot  that  operates  for  all 
aircraft  models. 

•  Model  the  disturbances  due  to  wind  gusts  and  refueling. 

•  Extend  QFT  design  approach  for  disturbances  at  the  output. 

•  Design  the  AFCS  using  QFT. 

•  Simulate  the  design  on  MATRIX,  and  EASY5x  to  validate  the  AFCS  design. 

1.7  Overview  of  the  Thesis 

This  thesis  describes  the  research  performed  in  the  development  of  the  AFCS. 
Chapter  II  discusses  QFT  and  the  expansion  of  the  QFT  equations  to  include  external 
output  disturbance.  The  air-to-air  refueling  AFCS  concept  is  discussed  in  Chapter  III. 
Chapter  IV  presents  the  aircraft  and  disturbance  models.  The  AFCS  design  process  is 
shown  in  Chapter  V'  followed  by  linear  and  nonlinear  simulations  in  Chapter  VI.  Finally, 
Chapter  VII  presents  the  thesis  conclusions  and  recommendations. 

1.8  Summary 

Future  USAF  cargo/transport  requirements  may  very  well  involve  long  flights  with 
many  air-to-air  refuelings  An  AFCS  to  automatically  pilot  the  aircraft  during  air-to-air 
refueling  will  greatly  reduce  the  pilot  workload  and  reduce  fatigue.  MIMO  QFT  design 
techniques  are  ideally  suited  for  the  design  of  such  a  flight  control  system  with  its  wide 
range  of  parameter  uncertainty.  This  thesis  lakes  on  the  challenge  of  validating  MIMO 
QFT  for  design  of  regulator  AFCSs.  And  demonstrating  automatic  air-to-air  refueling 
flight  control  systems  are  technically  achievable.  The  next  chapter  briefly  reviews  the 
QFT  process  and  the  development  of  the  output  disturbance  QFT  design  technique. 
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//.  QFT  and  Output  Disturbance  Rejection 


2.1  Introduction 

The  essence  of  QFT  designed  compensators  is  the  ability  to  accommodate  plant 

parameter  structured  uncertainty  and  reject  disturbances  [2],  This  Chapter  outlines  the 

QFT  method  and  how  nonlinear  plants  can  be  modeled  as  linear-time-invariant  (LTI) 

V 

plants  with  uncertainty.  First  a  brief  description  of  QFT  is  presented,  followed  by  the 
development  of  the  QFT  output  disturbance  model. 

2.2  Overview  of  QFT 

Quantitative  Feedback  Theory  is  a  powerful  design  method  for  synthesizing  fixed 
parameter  compensation  capable  of  controlling  a  system  in  the  realities  of  large  plant 
uncertainties,  disturbance  (rejection),  nonlinear  plant  models,  time  variations  of  plant 
parameters,  and  stringent  performance  specifications.  Included  in  QFT  is  the  ability  to 
design  a  compensator  to  cause  the  system  to  produce  outputs  that  meet  required  specifica¬ 
tions  based  on  particular  inputs  and  disturbance  rejection  requirements.  The  general  QFT 
design  process  is  based  on  the  diagram  presented  in  Figure  2.1,  which  represents  either 
a  multiple-input  single-output  (MISO)  or  a  MEMO  system  [2]  [3].  U  represents  the 
commanded  input  and  Dv  D2  are  disturbance  inputs.  G(s)  is  a  cascade  compensator  and 
F(s)  is  a  prefilter.  In  QFT,  (P(s)  represents  a  set  of  J  LTI  transfer  functions,  representing 
the  plant  at  different  operating  conditions.  The  designer  chooses  the  individual  transfer 
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Di  D2 


Figure  2.1  QFT  Controller  Design 


functions,  /^(s)  6  (P(s),  where  j=l,2,...,J,  so  that  the  entire  operational  envelope  or  the 
structured  uncertainty  of  the  plant  is  adequately  covered.  The  prefilter,  F(s),  is  designed 
so  that  the  system  meets  the  tracking  requirements  of  the  output  Y  to  the  input  U.  For 
a  regulator  control  system  (as  opposed  to  tracking)  the  commanded  input  is  zero,  thus,  no 
prefilter  is  required. 

The  objective  in  designing  the  G(s)  compensator  is  to  ensure  that  the  set  of  control 
ratios  (when  applicable): 


T  (s)  -  Y{s)  -  F(s)G(s)P(s) 
R  U^Cs)  1  ♦G(s)P(s) 


(2.1) 


and  the  disturbance  ratios: 


T  (S)  -  =  P(8) 

D  D,(s)  1  +  G(s)P(s) 


when  D2  =  0 


(2.2) 


T  fe)  -  Y(»)  ,  PCs) 

D  D2(s)  1  +  G(s)P(s) 


when  D,  =  0 


meet  the  desired  performance  and  disturbance  rejection  specification  for  all  Pj(s)  €  (P(s). 
The  nominal  loop  transmission  is  L0,  expressed  as 


) 
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L  =  GP 


(2.4) 


where  P„  represents  the  nominal  plant.  The  designer  selects  Pc  €  (P(s)  based  on  his/her 
preference  on  the  placement  of  the  templates'  boundaries.  To  form  the  stability  bounds, 
templates  are  derived  based  on  the  IP  information.  The  templates  are  plots  of  the  dB 
magnitude  versus  phase,  for  a  fixed  frequency  to,  for  each  P-}  contained  within  T  =  {Pj}. 
Lc  is  shaped  to  meet  the  tracking  BR(s)  and/or  disturbance  rejection  BD(s)  bounds  derived 
from  the  templates  on  the  Nichols  chart.  Combining  BR(s)  and  BD(s)  yields  the  optimal 
boundary  B0(jo)),  which  along  with  the  stability  contour  (U -contour),  constrain  Lc  such 
that  for  each  frequency  (Oj,  the  Lm  L^gi.)  point  must  lie  on  or  above  B0(jcOj),  and  not 
intersect  the  (/-contour.  Once  L„  is  formed,  G(s)  is  directly  derived  from  Eq.  (2.4). 

After  G(s)  is  formed,  whenever  applicable,  the  prefilter  F(s)  is  shaped  to  guarantee 
that  the  closed  loop  frequency  response  lies  withir  the  desired  tracking  specifications. 

In  the  case  where  -  0  (regulator),  the  prefilter  is  not  required  and  the  problem 
becomes  a  pure  disturbance  rejection  problem  .  In  this  thesis,  the  rejection  of  the  input 
disturbance  D2  at  the  system  output  is  the  desired  performance  specification. 

2.3  MIMO  QFT 

Multiple-input  multiple-output  QFT  design  method  requires  transforming  the  system 
into  equivalent  multiple-input  single-output  (MISO)  loops  as  shown  in  Figure  2.2  for  the 
3x3  case.  The  Schauder  fixed  point  Theorem  justifies  the  sequential  design  process, 
where  each  loop  is  individually  "closed"  [2],  Equation  (2.5)  represents  the  control  ratio 
.'y  which  relates  the  i*  output  to  the  f1  input  of  the  i,j  MISO  loop,  where  i  =  1,2,. ...m,  and 
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Figure  2.2  3x3  MISO  Equivalent  Loops 


j  =  l,2,...,m,/y  is  the  prefilter  of  the  i.j  MISO  loop,  /,  is  the  i*  loop  transmission,  and  d y 
represents  the  disturbance  input  in  the  i.j  loop.  The  are  the  reciprocal  of  the  terms 
contained  in  the  inverse  of  the  square  plant  matrix  P  =  [pti].  That  is,  P"  =  [p*J  =  P  \ 
where  [p‘y]  =  [1/<7U],  and  Q  =  [^].  The  disturbance  dv  given  by  Eq.  (2.6)  represents  the 
cross-coupling  between  loops.  The  disturbance  bounds  BD,  are  calculated  to  satisfy  Eq 
(2.7). 


f  1 .  +  da 

J  ij  i  u 

r  +  / ' 

i 


(2.5) 


",  -  -E 


it  =  1,  2 . m 


(2.6) 
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A  -  1,  2,  ...  m 


(2.7) 


1  W  U 


*  -  <7.* 


Note,  external  disturbances  are  not  accounted  for  in  this  development.  Only  cross- 
coupling  loop  to  loop  interaction  type  disturbances  are  considered.  A  separate  MIMO 
QFT  development  including  external  disturbances  is  presented  in  the  following  section. 

According  to  the  original  QFT  method  each  loop  is  designed  independently,  but  this 
leads  to  extreme  overdesign.  A  better  approach,  called  the  improved  method,  is  to  utilize 
the  known  structure  of  g ,  and  /s  of  the  loop  that  is  designed  first.  The  improved  method 
uses  the  known  structure  of  the  designed  loop  to  provide  a  more  accurate  estimate  of  the 
cross-coupled  disturbance  from  the  designed  loop  to  the  other  uncompensated  loops. 
Based  on  the  new  estimate  of  the  cross-coupled  disturbance,  new  disturbance  rejection 
bounds  for  the  uncompensated  loops  are  calculated.  The  improved  method  is  applied  after 
each  loop  is  compensated.  The  normal  QFT  loop  shaping  process  is  performed  for  a 
given  uncompensated  loop  using  the  new  disturbance  bounds.  The  improved  method 
reduces  the  overdesign  of  the  original  method  by  taking  into  account  that  compensated 
loop  will  have  less  cross -coupling  disturbance  intc  the  uncompensated  loops.  [3] 


2.4  MIMO  QFT  with  External  Output  Disturbance 

Output  disturbance  rejection  is  the  primary  design  criterion  in  this  thesis.  The 
previous  discussion  of  MIMO  QFT  did  not  consider  external  disturbance  in  the  calculation 
of  disturbance  rejection  bounds.  Therefore,  equations  including  the  external  disturbance 
are  developed  in  this  thesis.  The  following  discussion  outlines  this  development. 
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The  following  development  quantifies  external  uncertain  disturbances.  Figure  2.3 
represents  an  mxm  MEMO  closed-loop  system  in  which  F(s),  G(s),  P(s),  and  P,,(s)  are 
mxm  matrices.  (P(s)  =  { P(s) }  andlPd(s)  =  {Pd(s)}  are  sets  of  matrices  due  to  plant  and 
disturbance  uncertainties  respectively.  The  objective  is  to  find  a  suitable  mapping  that 
permits  the  analysis  and  synthesis  of  a  MIMO  control  system  by  a  set  of  equivalent  MISO 
control  systems. 


Figure  2.3  QFT  Controller  with  Output  External  Disturbance 
From  Figure  2.3,  the  following  equations  are  written 

y  =  P(s)u  +  P„(s)dCTt  u  =  G(s)v  v  =  r  -  y  (2.8) 

For  the  regulator  case  with  zero  tracking  input 

r  =  [o,  0,  o]T  (2-9) 
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From  Eqs.  (2.8)  and  (2.9)  where  henceforth  the  (s)  is  dropped  in  the  continuing  develop¬ 


ment 


v  =  -  y  u  =  -  Gy 


(2.10) 


which  yields 


v  =  -PGv  +  P ,d  .  (2.11) 

*  *  d  e*t 

Equation  (2.11)  is  rearranged  to  yield 

v  =[l  +  PG]  *P  d  ,  <212) 

Based  upon  unit  impulse  disturbance  inputs  for  de„,  the  system  control  ratio  relating  d^, 
to  y  is 


T,  =(l  +PG]'p<1  (2J3> 

Premultiply  Eq.  (2.13)  by  [I  +  PG]  yields 


[i  +PGJT,  -Pd  (214> 

when  P  is  nonsingular.  Premultiplying  both  sides  of  Eq.  (2.14)  by  P  1  results  in 


[P-'  +G]Td  =  P  'P,, 
Let 


(2.15) 


(2.16) 


2-7 


The  m:  effective  plant  transfer  functions  are  formed  as 


1  det  P 


the  Q  matrix  is  then  formed  as 


(2.17) 


</,« 

"  <7  |n 

1  Ip'u 

UP\  2  • 

■  1/P',. 

<?22 

-  7» 

~ 

l/p*2l 

l//»  • 

•  1/P  *2. 

7  mi 

<7  m2 

..  ii 

1  mm 

.  1  ip\, 

1/P'«2  ’ 

■  i/pV 

where  P  =  [/;„],  P  1  =  [/>'„]  =  (1  /<?,,].  and  Q  =  [qtJ]  =  [ l//>“y)  Partition  the  P '  matrix  as 
follows; 

P"  =[/>;  ]  =[  1/7,  1  =A  +B  <219> 

where  A  is  the  diagonal  part  of  P 1  and  B  is  the  balance  of  P Thus  Xy  =  \/qu  =  p\,  bu 
=  0,  and  b{j  =  \/q0  =  p\s  for  i  *  j.  Substituting  Eq.  (2.19)  into  Eq.  (2.15)  with  G  diagonal, 
results  in 

[  A  +  B  +  G  ]Td  =  [  A  +  B  ]Pd  (2-2°) 

Rearranging  Eq.  (2.20)  produces 

Td  =  [  A  +  G  ]' '[  APd  +  BPd  -  BTd  ]  (2  21) 

'This  equation  defines  the  desired  fixed  point  mapping,  where  each  of  the  nr  matrix 
elements  on  the  right  side  of  Eq.  (2.21)  are  interpreted  as  MISO  problems.  Proof  of  the 
fact  that  the  design  of  each  MISO  system  yields  a  satisfactory  MIMO  design  is  based  on 
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Figure  2.4  3x3  MISO  Equivalent  Loops  for  External  Output  Disturbance 

the  Schauder  fixed  point  theorem.  [2]  The  theorem  defines  a  mapping  Y(Td) 

Y(Td)  *  [  A  +  G  V[  APd  ♦  BPd  -  BTd  ]  (2.22) 

where  each  member  of  Td  is  from  the  acceptable  set  3d.  If  this  mapping  has  a  fixed  point, 
i.e.,  Td  €  3d,  then  this  Td  is  a  solution  of  Eq.  (2.21). 

Figure  2.4  shows  the  effective  MISO  loops  resulting  from  a  3x3  system.  Since  A  and 
G  in  Eq.  (2.21)  are  diagonal,  the  (1,1)  element  on  the  right  side  of  Eq.  (2.22)  for  the  3x3 
case,  for  a  unit  impulse  input,  provides  the  output 
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Sr 

II 

Pa  Pd  P, 

( 

jq  (2.23) 

1 

^11  $12  ^13 

,  ^12  ^13  , 

Equation  (2.23)  corresponds  precisely  to  the  fust  structure  in  Figure  2.4.  Similarly, 
each  of  the  nine  structures  in  Figure  2.4  corresponds  to  one  of  the  elements  of  Y(Td)  of 
Eq.  (2.22).  The  control  ratios  for  the  external  disturbance  inputs  d a  and  the  correspond¬ 
ing  outputs  y j  for  each  feedback  loop  of  Eq.  (2.22)  have  the  form 

v  -  k '(d)  (2.24) 

* 

where  w,;  -  qj(l  +  gqO  and 


x  P  m  1 

de  =  E  —  -  E  - 

•  * q*  kt,  q* 


x  =  number  of  disturbance  inputs 
m  -  dimension  of  square  MIMO  system 


(2.25) 


Thus  Eq  (2.25)  not  only  contains  the  loop  interaction  but  also  the  external  disturbances. 

Additional  equations,  quantifying  the  internal  cross -coupling  due  to  external 
disturbances,  are  derived  to  utilize  the  improved  method  QFT  design  technique.  These 
equations  are  used  to  define  the  disturbance  bounds  for  subsequent  loops  based  on  the 
completed  design  of  a  single  loop.  For  this  development,  the  equations  for  the  case  of 
a  2x2  MIMO  system  are  presented. 

From  Eq.  (2.24)  and  for  the  1-2  loop  case,  which  is  the  output  of  loop  1  due  to 
disturbance  input  2,  including  the  cross -coupling  terms  from  loop  2,  yields,  for  unit 
impulse  inputs,  the  following  control  ratio 
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V,  =y12  =wi,(  dJ  = 


^11 


1  +  L, 


Prt  P„ 

dH  +  <*33 

9n  ^12 


S. 

Ql2 


(2.26) 


Substituting  in  for  t 


V  - 


1  +  L, 


P.  P.  q,,d 

dL)  +  ^  »a 

9)1  Q)2  (i  +  ^2)^12 


(2.27) 


t.  = 


9n 


1  +  L, 


+  ^2^12  +  ^(^  +^2^1i  ^22^1  lde 

(1  +  L2)quq12 


(2.28) 


rearranging,  and  substituting  in  for  dc^  yields 


(1  +  ^(^,^12  +^>dD<lll^  “  ^22^1 


V,  ‘  “ 


pd  pd. 

d.i  + _ 


1% 


122 


_*d 

q2.j 


(2.29) 


(1  i-L.Kl  +  L2)qJ2 


(1  +L2)(Pduq,2  -rP^q,,) 


V.  = 


^22^,.^,, 

q2. 


-%Pa 


^22^11^ 

q21 


(2.30) 


(1  +  L,)(l  +  L2)q12 


(1  +  ^2)(Pd, ,9|2  +  P^ll^l  ~  tl22<qilPd[J  ~  ^11^21^  +  922^,1^ 

=  "(T+L.Xi  +  L2)q12q2i 


(2.31) 


where  y12  =  q!  ^2^12^12-  Solving  for  t^  yields 
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((1  +  L  )(1  +  L,))  = 


(1  +  L,)(P  q12  +  P,qn) 


^ll^d 


in 


y  n^«j„ +  y  n^d 


m2 


(3.32) 


t„«l  -  L,)d  +L2)  -Yn)  = 

(i  *<i  .ljp^  -rlAu 

Hl2  4l2 


(3.33) 


ld..  = 


L,-qJip  .(l  *L,  -  y „)Pj 

q.2 _ 

L,(  1  +  L2)  +1  +  L2  -y12 


(3.34) 


Equation  (3.34)  is  rearranged  as  follows: 


112  +  p 

.  1  +L2  ~  Y  12 _ * 

j  „  L,q  *W 

1  +L2  "Yn 


(3.35) 


From  Eq.  (3.35)  the  effective  plant  is  defined  as 


q 


u. 


qua  + 

1  +L2  'Yn 


(3.36) 


Substituting  Eq.  (3.36)  into  Eq.  (3.35),  yields 

Thus,  in  general,  for  the  2x2  case,  the  improved  method  control  ratio  of  the  f1 
external  disturbance  input  to  the  /th  system  output  is 
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(3.37) 


t 


q„a 

i  +g,q 


+  p. 


12 


'*  = 


9.d  +^*) 


1  +^u 


where  /  -1,2  and  k  *  i 


(3.38) 


The  new  disturbance  bounds  are  calculated  at  a  given  frequency  to  satisfy 


Bd  * 


LJjX 

9.0  +  ^) 


1  +  «,9. 


where  /  =  1,  2  and  k  *  i 


(3.39) 


or 


1  *W.J27 


fA‘i. 

9.d  +  ^) 


where  i  -  1,  2  and  /;  *  /  (3.40) 


The  improved  QFT  method  uses  these  equations  to  reduce  overdesign  inherent  in  the 
original  design  process.  The  order  in  which  loops  ate  designed  is  important.  Any  order 
can  be  used,  but  some  orders  produce  less  overdesign  (less  bandwidth)  than  others.  The 
last  loop  designed  has  the  least  amount  of  overdesign,  therefore  the  most  constrained  loop 
is  done  first  by  the  original  method.  Then  the  design  is  continued  through  the  remaining 
loop  by  the  second  method.  The  first  loop  is  then  redesigned  using  the  improved  method. 

At  this  point  it  is  important  to  point  out  that  when  the  disturbance  rejection 
specification  is  considered,  the  designer  must  decide  how  much  disturbance  rejection  is 
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for  cross-coupling  and  how  much  is  for  external  disturbances.  In  other  words  the  designer 
can  "tune"  the  disturbance  rejection  specification  depending  on  the  nature  of  the 
disturbance  for  a  particular  loop.  For  example,  if  one  loop  is  only  effected  by  external 
distrubance,  the  disturbance  rejection  specification  would  consider  external  disturbance 
only.  But  if  the  loop  disturbance  was  a  mix  of  cross-coupling  and  external  disturbance, 
the  designer  would  have  to  "tune"  the  disturbance  rejection  specification  accordingly. 
Since  each  loop  may  not  exhibit  the  same  disturbance  characterisitics,  disturbance 
rejection  specification  tuning  provides  flexibilty  in  the  QFT  design  process  . 


2.5  Summary 

This  chapter  presented  the  basic  improved  QFT  method  and  its  extension  to  MIMO 
systems.  The  development  to  include  the  effects  of  external  output  disturbances  is 
detailed.  Equations  to  calculate  the  new  disturbance  bounds  for  the  improved  method  are 
discussed.  The  next  chapter  covers  modeling  the  aircraft  and  the  external  disturbances. 
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III.  Air-to-Air  Refueling  FCS  Design  Concept 


3.1  Introduction 

The  aircraft  modeled  in  this  thesis  is  the  cargo  variant  of  the  135  class  aircraft,  the 
C-135B.  This  aircraft  is  chosen  because  of  the  availability  of  the  aerodynamic  data  [1], 
This  chapter  describes  the  modeling  of  the  C-135B.  A  Mach-hold,  altitude  hold,  and 
wing-leveler  autopilot  is  included  in  the  C-135B  model  [6].  Wind  gusts  and  fuel  transfer 
disturbance  models  are  developed  in  this  chapter,  as  well  as  the  AFCS  concept. 

3.2  C-135B  Modeling 

EAS  Y5x  is  used  to  develop  the  state-space  six-degrees-of-freedom  bare  (uncontrolled) 
aircraft  model.  EASY5x  is  a  computer  aided  design  (CAD)  tool,  written  by  Boeing 
Computer  Services,  used  to  model,  simulate,  and  analyze  dynamic  systems  [5]  EASY5x's 
graphical  interface  expedites  modeling  and  analysis.  Generic  pre- written  air  vehicle 
modules  further  simplify  modeling  in  EASY5x.  The  user  need  only  provide  aircraft 
stability  derivatives,  flight  conditions,  and  desired  input/outputs. 

The  basic  aircraft  model  is  illustrated  in  Figure  3.1.  The  LO  block  calculates  the 
longitudinal  forces  and  moments,  while  LD  does  the  same  for  the  lateral-directional  axes. 
The  AV  block  includes  the  aerodynamic  variables  and  calculates  changes  in  speed  and 
rates  due  to  wind  gusts.  Also,  it  computes  dynamic  pressure  and  angle-of-attack  (AOA). 
The  SD  block  contains  the  six-degrees-of-freedom  equations  of  motion.  It  calculates 
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Figure  3.1.  Bare  Aircraft  Plant 


attitude,  body  axes  rates,  Euler  angles,  and  body  axes  angular  rates  based  on  data  supplied 
by  the  LO  and  LD  blocks.  Further  information  can  be  obtained  from  the  EASY5x  user 
manual  [5]. 

Sixteen  bare  aircraft  plants  are  developed  to  account  for  the  uncertainty  of  the  C- 
135B  during  air-to-air  refueling.  The  16  models  are  based  on  two  different  coefficients 
of  lift,  CL  =  0.2,  0.6,  for  eight  different  aircraft  weights,  each  flying  Mach  0.69  at  28,500 
feet.  These  discrete  values  are  selected  based  on  the  availability  of  data,  normal  refueling 
speed  and  altitude,  and  represent  weights  ranging  from  empty/low  fuel  to  loaded/full  fuel 
aircraft.  Typically,  during  air-to-air  refueling,  the  C-135B  will  have  a  CL  between  0.27 
and  0.45  [1].  Therefore,  the  16  plant  models  envelop  the  structured  uncertainty  of  the  C- 
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135B  during  air-to-air  refueling.  Tables  containing  the  aerodynamic  stability  derivatives 
for  all  16  plant  models  are  provided  in  Appendix  A. 

The  six-degrees-of-freedom  state-space  models,  generated  by  EASY5x,  are  loaded 
into  MATRIX,.  MATRIX,  is  a  programmable,  matrix  calculator  with  graphics. 
MATRIX,  was  designed  by  control  engineers  to  aid  in  the  design  of  complex  control 
problems.  It  is  ideally  suited  to  manipulate  matrices  and  apply  both  input/output 
(classical)  control  and  state-space  (modem)  control  [7], 

MATRIX,  is  used  to  design  the  autopilot  using  root-locus  design  techniques.  The 
autopilot  is  designed  to  control  all  16  plant  cases.  The  bare  aircraft  and  autopilot  are 
shown  in  Figure  3.2.  An  autopilot  is  used  for  two  reasons:  (1)  autopilots  reduce  the  high 
frequency  cutoff  of  the  aircraft  (2)  all  aircraft  have  autopilots.  Lowering  the  cutoff  of  the 
aircraft  reduces  high  frequency  parameter  uncertainty.  Since  autopilots  are  available, 
using  it  in  the  QFT  design  eliminates  duplication  of  a  control  system  to  provide  input  to 
the  bare  aircraft,  reducing  cost  and  overhead.  The  inputs  to  the  autopilot  are  thrust, 
elevator,  aileron,  and  rudder  commands.  The  outputs  are  z-position  (altitude),  x-position, 
and  y-position  in  a  local  inertial  reference  frame  where  x  is  positive  out  the  nose  of  the 
aircraft,  y  out  the  right  wing,  and  altitude  is  positive  up.  The  three  outputs  frame  of 
reference  is  translated  from  the  aircraft  center  of  gravity  (eg)  to  the  approximate  location 
of  the  air-to-air  refueling  receptacle  on  the  top  of  the  aircraft. 

In  order  to  have  a  square  plant,  3x3,  one  of  the  autopilot  inputs  is  not  used.  The 
Mach  hold  command  input  is  used  to  control  the  x  position,  altitude  hold  controls  altitude, 
and  the  rudder  command  is  used  to  control  the  y  position.  Mach  and  altitude  are  self 
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Figure  3.2  C-135B  Bare  Aircraft  with  Autopilot 


evident,  rudder  is  chosen  over  aileron  because  the  rudder  does  not  roll  the  aircraft  as 
much  as  the  ailerons.  By  using  the  rudder  for  the  QFT  controller,  good  performance  is 
obtained  while  leaving  the  aileron  controller  to  handle  wing  leveling.  See  Appendix  C 
for  aircraft  response  to  step  input  commands  on  all  four  autopilot  command  inputs. 

3.3  Disturbance  Modeling 

Disturbance  models  are  generated  by  developing  augmented  state-space  models  of  the 
aircraft  in  the  presence  of  wind  gusts  and  fuel  transfer  inputs.  This  development  starts 
by  examining  the  aerodynamic  moments  and  forces. 

3.3.1  Pitch  Plane  Wind  Induced  Disturbance 

In  the  pitch  plane,  the  total  aircraft  airspeed  is  a  combination  of  wind  velocity 
and  aircraft  velocity  in  the  x  and  z  stability  axis  directions. 
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Airspeed  - 


(3.1) 


UA 

WA 


0  +  u  +  uw 
w  +  w 

w 


Where  U  is  the  steady-state  aircraft  speed,  u  &  w  are  the  perturbation  speed  in  the 

x  and  z  directions  respectively,  and  u*  &  w».  are  the  aircraft  speeds  due  to  wind  in  the 
respective  x  and  z  directions. 

The  pitch  plane  aerodynamic  moments  and  forces  are  in  the  form  of  the  following 
functions 


where 


u  qSc  r  .  ,..,180. 

M  =  Cmra.  q ,  5)( - ) 

71 

F  =  ®  C  (a,  6) 
m 


c  As  rx,180. 

F  =  -i—  C(a,  q,  6)( - ) 

mU  n 


*WA2) 


a 


=  arctan 


\ 


(3.2) 


(3.3) 


and  the  a  effects  in  the  M  and  Fz  equations  are  neglected. 

It  is  important  to  note  that  the  following  equations  of  motion  are  written  with  respect 
to  the  inertial  kinematic  variables  q,  U,  and  W.  However,  the  perturbations  in  U  and  W 
are  u  +  uw  and  w  +  w„  respectively.  The  q  perturbation  is  q  +  qw. 
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The  q  equation  of  motion  is 


q  =  Mao.  +  Mqq  +  Maw  +  Mhb 


t  iEELc.[UA(u  ♦  u.)  ♦  WA(w  *  w,.)j  *  w. 


=  M  a  +  Mo  +  M  q  +  M.b 

a  q  1  q 1  h*  0 


+  _180_  P5£_CJ  (0  +  u  +  u J(u  +  uw)  +  (w  +  wj2  +  Ma~ v 
n  I,  ^  0 


assuming  u,  w,  uw,  and  ww  «  0  Eq.  (3.4)  reduces  to 


q  =  Maa  +  Mqq  *  Mqqw  +  Ma6  +  Mu[  u  +  u  J  +  Afa— wv 


grouping  the  wind  terms  results  in 


q  =  Ma  a  +  Mq  +  Mbb+  Mu  *  [  M  a  +  A/uw  +  M— ww ) 


A  similar  derivation  for  the  a  and  u  equations  of  motion  yield 


1  tfw.  1 

a  =  Za  +  Zo  +  Z  u  +  Z.6  +  [Z  u  +  (Z  -  l)o  + - +  — Z  w  ] 

a  y*  w  6  1  u  w  x  y  'T*  q  ^  ^  a  wJ 


u  *  *aa  ♦  Xqq  *  Xuu  *  Xbb  +  [Zw,  ♦  Xqqw  +  ±X*  J 


Hence,  the  effect  of  a  moving  airmass  is  akin  to  "process  noise",  i.e.,  if  the  state  is 


x  -  (0,  u.  a,  q) 
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the  state-space  equation 


x  =  A  x  +  B  u 


is  augmented  to  yield 


(3.9) 


(3.10) 


where 
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(3.11) 


and  the  disturbance  vector  is 


^  pilch 


Uw 

w 


(3.12) 


\  dw 

Note,  - ,  which  is  an  angle-of-attack  rate,  is  neglected.  This  simplification  is  in 

0  di 


line  with  discarding  of  the  &  effects.  Further  simplification  is  achieved  by  neglecting 
<7W  .  Hence,  rpiut  and  dpilch  in  the  pitch  channel  are  reduced  to 
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J.5.2  Lateral  Channel  Wind  Induced  Disturbance 

For  the  lateral  channel,  arguments  similar  to  the  pitch  channel  development  are 
made  to  generate  the  augmented  state-space  equation  for  the  effects  of  wind.  In 
conformity  with  the  level  of  modeling  in  the  pitch  channel,  if  the  state  is 

x  =  (y,  p,  r,  p )  (314) 

rta  and  d^  are  derived  to  be 


(315) 


3.3.3  Disturbance  Due  to  Refueling 

As  the  aircraft  is  taking  on  fuel,  the  moments  of  inertia  and  eg  are  changing. 
These  variations  alter  the  flying  behavior  of  the  aircraft.  This  change  in  aircraft  behavior 
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is  modeled  as  a  disturbance.  Based  on  first  order  approximation,  refueling  (rf)  caused 
disturbance  is  confined  to  the  pitch  channel.  Again,  paralleling  the  pitch  channel 
development,  the  state-space  model  is  augmented.  rri  and  drf  are  found  to  be 


0 

0 

12  32.2 

- *  c* - 

100  I„ 

J_* 32  2 

0  m 


dK  =  [fuel  transfer  rate  (lbs/min)] 


(3.16) 


Where  12/100  is  the  percentage  of  the  mean  aerodynamic  cord  the  eg  is  allowed  to  shift 
during  refueling  and  m  is  the  mass  of  the  aircraft. 

3.3.4  Total  Disturbance 

The  total  modeled  disturbance  is 

rd  =  r  „  hd  „  +  r  d^,  +  r\drf  (3.i7) 

pilch  pilch  Ut  Ul  rf  rf 

where  rpil£h  and  rrf  augment  the  state-space  equation  containing  the  states  that  identify 
pitch  plane  flight  behavior.  In  the  same  manner,  T,*  augments  the  lateral  channel  states. 
The  result  of  the  disturbance  summation  for  the  augmented  disturbed  states  (0,  u,  a,  q, 
p.  r,  p)  is  shown  in  Eq.  (3.18). 
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The  state-space  equation  now  takes  on  the  form 

x  =  Ax  +  Bu  +  Td  (3.19) 

y  =  C  x 

Note  as  the  disturbances  are  applied  to  the  bare  aircraft,  they  enter  into  the  inner 
most  loop,  around  which  the  autopilots,  and  later,  the  QFT  loops  are  closed. 


3.4  Plant  and  Disturbance  Matrices 

Based  on  zero  initial  conditions,  then  from  Eqs  (3.19) 

sx  =  Ax+Bu  +  rd  (3.20) 

»-[.l-A]-,B.4l-A]"rd  <32‘) 


3-10 


y  =Cx  =  c[s!  -  a]  '  Bu  +  c[sl  -  A  ]  '  Td 


(3.22) 


where 


=  P(s)u  +  Pd(s)d 


P(s)  =  c[sl  -  A  ]  '  B 
Pd(s)=c[sl-A]‘‘r 


(3.23) 


Equation  (3.22)  is  represented  in  Figure  2.3.  Thus  the  QFT  formulation  of  Section 
2.4  is  applicable  for  this  problem. 


3.5  Control  Problem  Approach 

The  tanker’s  position  ia  assumed  fixed  and  hence  the  receiver  aircraft's  position  is 
measured  from  this  frame  of  reference.  In  this  approach  the  control  problem  can  be 
viewed  as  a  formation  flying  problem.  The  receiver  maintains  the  total  obligation  of 
regulating  its  position.  The  tanker  is  free  to  change  course,  altitude,  and  velocity  while 
the  receiver  compensates  for  these  changes  and  maintains  position. 

Equations  are  developed  that  identify  perturbations  from  the  set  position.  These 
perturbations  are  viewed  as  disturbances  by  the  receiver.  The  perturbations  are  caused 
by  wind  gusts  and  disturbances  due  to  refueling.  Other,  unmodeled,  disturbances  may 
include  the  tanker  changing  course.  The  control  problem’s  goal  is  to  minimize  the 
perturbations  to  be  within  a  specified  volume  of  space  where  the  refueling  boom  can 
operate.  Normal  boom  operating  position  and  length  defines  this  volume. 
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The  KC-135  tanker  refueling  boom  has  the  following  operational  constraints,  (1) 
nominal  boom  operation  position  is  30  degrees  down  from  horizontal,  (2)  the  boom  can 
move  as  much  as  four  degrees  up  and  down  from  normal  position  and  continue  delivering 
fuel,  (3)  it  can  move  as  much  as  ten  degrees  up  and  down  from  normal  position  and 
maintain  its  connection  to  the  receiver,  but  cannot  deliver  fuel,  (4)  horizontal  movement 
is  limited  to  10  degrees  left  and  right  while  maintaining  fuel  flow,  (5)  the  disconnect  limit 
horizontally  is  15  degrees  left  and  right,  (6)  nominal  boom  length  is  477.5  inches  (39.8 
ft),  (7)  it  can  expand  or  constrict  13.5  inches  and  maintain  refueling,  (8)  it  can  expand 
or  constrict  as  much  as  73.5  inches  and  maintain  contact  but  not  refueling.  [1] 

These  dimensions  provide  a  maximum  perturbation  from  nominal  boom  position  of 
approximately  2.85  feet  up  or  down,  7  feet  left  or  right,  in  order  to  maintain  fuel  flow. 
In  order  to  maintain  connection,  the  maximum  perturbation  can  be  7.5  feet  up  or  down, 
and  11.5  feet  left  or  right. 

Using  the  tanker  as  the  point  of  reference,  the  relationship  in  Figure  3.3  can  be  used 
to  develop  the  equations  required  to  define  the  regulation  control  problem.  R  is  the 
nominal  boom  length  measure  from  the  boom  hinge  point  on  the  tanker.  Z  is  the  vertical 
distance  between  the  boom  hinge  point  and  receiver  aircraft's  refueling  receptacle.  X  is 
the  horizontal  distance  between  these  same  points.  Y  measures  the  distance  between  the 
center  line  of  the  boom  hinge  point  and  receiver  receptacle. 

The  following  equations  arc  derived 

R2  =  X2  +  Y2  +  Z2  (3.24) 
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where 


R  =  R  +  r 

x  =  *  +  *  (3.25) 

Y=  Y  +  y 

Z  =  Z  +  z 

which  are  the  sums  of  the  nominal  positions  (overbar  terms)  and  perturbations  (lower  case 
terms).  Substituting  Eq.  (3.25)  into  Eq.  (3.24)  and  squaring  the  terms,  yields 

R2  +  2rR  +  r 2  =  X2  +  2xX  +  x2  +  Y2  +  2y?  +  y2  +  22  +  2z2  +  z2  <3  26) 


Figure  3.3  Control  Problem  Geometry 


Since  r,  x,  y,  z  «  R,  X,  2  respectively,  Eq.  (3.26)  is  approximated  as 

ft2  +  2rR  =  X2  +  2xX  +  ?2  +  2y?  +  22  +  2z2  <3-27) 

Taking  the  derivative  with  respect  to  time  where  the  overbar  terms  are  constant  yields 
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(3.28) 


x  =  Xt-xr 

y  =  yT-yR 


T  -  Tanker 
R  -  Receiver 


(3.30) 


Thus,  r  =  0  if  and  only  if,  x  =  y  =  z  =  0.  Therefore,  the  control  problem  is  to  design  the 
compensator,  G,  of  Figure  3.4  that  will  satisfy  Eqs.  (3.29)  and  (3.30). 


Figure  3.4.  Control  Problem 
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3.6  Summary 

The  C-135B  model  is  developed  using  the  CAD  program  EASY5x.  An  autopilot  is 
designed  for  the  bare  aircraft  model.  The  QFT  plants  are  the  16  aircraft  models  with 
autopilot.  External  disturbances  enter  the  plant,  the  aircraft  state-space  model  is 
augmented  to  model  the  disturbance.  The  regulation  control  problem  is  identified.  The 
next  chapter  is  devoted  to  the  generation  of  the  LTI  plants  and  disturbances. 
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IV.  LTI  Plant  and  Disturbance  Model  Generation 


4.1  Introduction 

Before  QFT  output  disturbance  rejection  is  applied  to  the  design  of  an  AFCS,  the  LTI 
plant  matrix,  P(s),  and  disturbance  matrix,  Pd(s)  must  be  generated.  This  chapter  presents 
the  results  of  applying  the  methods  of  chapters  11  and  111  in  forming  the  plant  q's.  First, 
the  state-space  models  are  generated  and  transformed  into  LTI  transfer  functions 
representing  the  16  plants.  Next,  the  q's  are  formed  by  applying  the  plant  inversion 
process.  Finally,  disturbance  transfer  functions  are  generated  for  all  16  plants.  The 
results  represent  the  uncertainty  in  the  plants  as  well  as  the  disturbance. 

4.2  Plant  Transfer  Function  Generation 

First  the  bare  aircraft  state-space  model  is  generated  in  EASY5x.  This  model  is 
loaded  into  MATRIX*  where  the  state-space  model  with  autopilot  is  formed.  The  plant 
matrix  P(s)  of  Eq.  (3.23)  is  calculated  in  the  MIMO/QFT  CAD  package  from  the  state- 
space  model. 

x  =  A  x  +  B  u  (4  i) 

y=Cx 

The  transfer  functions  are  calculated,  as  shown  in  Sec.  3.4,  from  the  following  relationship 
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y-C(sI-A)4Bu  (4  2) 

y(s)  =  P(s)u(s) 

Therefore 

P(s)  =  C(sI-  A)hB  <4-3) 

For  the  3x3  MIMO  system,  P(s),  is  a  matrix  of  nine  transfer  functions.  The  transfer 
functions  along  the  matrix  diagonal,  Pu(s),  P22(s),  and  P33(s),  represent  the  respective 
output  response  to  the  given  input,  e  g.,  output  1  to  input  1.  Off-diagonal  terms  represent 
the  cross-coupling  between  channels,  e.g.,  output  1  to  input  2.  The  columns  indicate 
input,  column  1/input  1,  and  the  rows  are  the  outputs,  row  2/output  2.  Hence,  the  MISO 
loop  {2,1}  is  the  output  channel  2  due  to  input  channel  1.  As  seen  in 
Figure  4.1,  Figure  4.2,  the  dominant  input  for  output  1  is  input  2,  see  MISO  loop  {1,2}. 
This  indicates  a  strong  correlation  between  channels  1  and  2. 

This  correlation  is  expected.  Channel  1  is  the  altitude  hold  autopilot  and  channel  2 
is  the  Mach  hold.  The  altitude  is  controlled  by  using  the  thrust.  To  increase  altitude  the 
throttle  is  advanced,  which  increases  the  speed  and  hence  lift,  therefore  increasing  the 
altitude.  The  reverse  is  true  for  decreasing  altitude.  The  speed  is  controlled  with  the 
elevator  that  controls  the  pitch  angle  of  the  aircraft.  To  increase  speed  the  aircraft  is 
pitched  down,  decreasing  the  drag,  resulting  in  an  increase  of  aircraft  airspeed.  The 
opposite  is  true  to  decrease  airspeed.  Lift  and  drag  are  inversely  related,  and  thus  highly 
coupled.  An  increase  in  drag  decreases  lift.  Therefore  a  decrease  in  airspeed  (increased 
drag),  also  results  in  a  decrease  in  altitude  (reduced  lift).  As  the  Mach  hold  autopilot 
controls  airspeed  it  adjusts  pitch  angle,  increasing  or  decreasing  drag.  This  results  in 
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Figure  4. 1  P(s)  Log  Magnitude  Plot 


changes  of  altitude,  which  is  controlled  by  the  thrust.  The  thrust  is  adjusted  using  the 
engines  which  have  a  longer  lag  time  (longer  reaction  time)  than  the  elevators.  This 
causes  a  delay  between  when  the  altitude  hold  senses  change  in  altitude  and  the  altitude 
correction.  Thus,  the  correlation  of  input  2  to  output  1. 
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Notice  next  that  channel  3  (y  position)  is  essentially  decoupled  from  channels  1  and 
2,  MISO  loops  {1.3}.  {2,3}.  {3.1},  and  {3,2}.  Again  this  is  expected.  Channel  3  is  the 
lateral  channel,  derived  separately  from  channels  1  and  2,  longitudinal  (pitch)  channel. 
Decoupling  channel  3  allows  the  3x3  MISO  system  to  be  broken  down  to  two  design 
problems,  a  2x2  MISO  system  and  a  lxl  SISO  system.  Thus  the  QFT  improved  method 
is  applied  only  to  channels  1  and  2. 

The  single  dotted  line  in  the  MISO  loops  of  Figure  4.1  is  the  log  magnitude  plot  of 
the  disturbance  rejection  specification.  Log  magnitude  of  P's  above  this  line  indicates  that 
the  output  due  to  cross-coupling  disturbance  exceeds  the  disturbance  rejection  specifica¬ 
tion  for  the  frequencies  that  lie  above  the  line.  Note  that  P(s)  MISO  loops  {1,3},  {2,3}, 
{3,1},  and  {3,2}  aie  below  the  disturbance  rejection  line,  reinforcing  replacement  of  the 
3x3  system  by  the  2x2  and  lxl  systems  as  discussed  previously. 

Finally  notice  in  Figure  4.1  the  tight  qrouping  of  the  16  plants  over  the  frequency 
range  plotted.  This  indicates  little  uncertainty  in  magnitude  for  the  16  plants.  This  is  an 
added  robustification  benefit  due  to  the  "inner  loop"  autopilots. 

The  QFT  Q(s)  matrix  is  formed  in  the  MIMO/QFT  CAD  program  by  matrix  inverting 
the  transfer  functions  of  the  P^s)  matrix.  The  q y  transfer  functions  for  the  16  plants  are 
shown  in  Appendix  B.  Figure  4.2  shows  the  log  magnitude  (Lm)  plots  of  the  q4  transfer 
functions  for  all  plant  cases.  As  mentioned  above,  channel  3  is  consider  completely 
decoupled  from  channels  1  and  2.  Thus,  for  design  purposes  ql3,  q^,  q3t,  and  q32  arc  set 
to  zero  since  their  cross-coupling  contribution  is  essentially  zero. 
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Figure  4.2  Q(s)  Log  Magnitude  Plots 


4.3  Disturbance  Transfer  Function  Generation 

The  disturbance  models  Pd(s)  of  Eq  (3.23)  are  calculated  similar  to  P(s).  The 
disturbance  augmented  state-space  model  is  developed,  as  detailed  in  Chapter  ///. 
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x  =  A  x  Bu  +  Td 
y  =  Cx 


(4.4) 


The  disturbance  matrix  Pd(s)  is  given  by 

Pd(s)  =  C(sl  -  A)'T  (4.5) 

Pd(s)’s  are  calculated  for  all  16  plants.  Like  P(s),  Pd(s)  is  a  matrix  of  transfer 
functions  with  the  same  matrix  diagonal,  off-diagonal  characteristics. 

The  disturbance  driving  functions  are:  wind  gusts  in  the  vertical  up  direction, 
horizontal  left  to  right,  and  fuel  flow  rate.  The  magnitude  of  the  wind  gusts,  horizontal 
and  vertical,  are  10  feet/sec,  applied  in  a  step  function.  The  fuel  flow  rate  is  10,000 
pounds  per  minute  (maximum  fuel  delivery  rate  of  the  KC-135  tanker),  applied  in  a  ramp 
function 

Figure  4.3  presents  the  Lm  Pd(s)  plots  for  all  plant  cases.  It  is  important  to  keep  in 
mind  that  the  wind  and  fuel  disturbances  enter  the  bare  aircraft  but  the  Pd(s)  models 
include  the  autopilots.  Therefore  the  Pd(s)  output  will  not  only  exhibit  the  effects  of  the 
disturbance  input  but  also  the  cross-coupling  nature  of  the  autopilot  discussed  previously. 

In  Figure  4.3  column  1  maps  disturbance  input  1  (refueling)  to  the  three  system 
outputs,  altitude,  x  position,  and  y  position  respectively.  Column  1  shows  very  little 
effect  from  refueling  on  the  outputs  (MISO  Loops  {1,1),  {2,1]  and  {3,1}).  Column  2 
demonstrates  that  disturbance  input  2  (vertical  wind)  has  a  strong  effect  on  both  altitude 
and  x  position.  Column  3,  input  3  (horizontal  wind)  effects  only  y  position. 

In  Figure  4.3  notice  that  MISO  loops  {*,1 }  and  {*,2}  show  two  distinct  groupings 
of  Pd(s)  log  magnitude  plots.  Each  grouping  represents  the  plant  differences  based  on  CL. 
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MISQ  loop  \  1.3  [ 


Figure  4.3  Pd(s)  Log  Magnitude  Plots 


The  distinct  groupings  indicate  that  the  external  disturbances  effect  the  CL  =  0.6  plants 
more  than  the  CL=  0.2  This  is  also  evident  in  the  Pd  time  response  plots  presented  in 
Appendix  D. 

The  single  dotted  line  in  the  MISO  loops  of  Figure  4.3  is  the  log  magnitude  plot  of 
the  disturbance  rejection  specification.  Log  magnitude  of  Pd's  above  this  line  indicates 
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that  the  output  due  to  disturbance  exceeds  the  disturbance  rejection  specification  for  the 
frequencies  that  lie  above  the  line.  Pd(s)  MISO  loops  {i,3}.  {2,3},  {3,1},  and  {3,2}  are 
below  the  disturbance  rejection  line,  thus,  further  justifying  the  break  down  of  this  3x3 
system  into  2  systems,  2x2  and  lxl. 

4.4  Summary 

In  this  chapter,  LTI  models  and  Bode  plots  for  P(s)  and  Pd(s)  are  generated.  The 
calculation  of  the  Q(s)  matrix  is  discussed  as  well  as  the  significance  of  channel  3  being 
uncoupled  from  channels  1  and  2.  The  next  chapter  describes  the  AFCS  QFT  design. 
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V.  QFT  AFCS  Design 


5.1  Introduction 

This  chapter  details  the  QFT  AFCS  design.  First  the  disturbance  rejection 
specification  is  identified.  Next,  design  of  the  loop  transmissions  for  all  three  channels 
is  described.  This  chapter,  nor  this  thesis,  is  a  detailed  step-by-step  guide  for  the  QFT 
design  process.  The  reader  is  assumed  to  have  at  least  a  basic  understanding  of  the  QFT 
method. 

5.2  Disturbance  Rejection  Specification 

The  primary  goal  in  designing  the  AFCS  system  is  to  regulate  the  position  of  the 
aircraft  receiving  fuel  relative  to  the  tanker.  As  discussed  in  Section  2.4,  any  deviation 
from  the  nominal  set  position  is  considered  a  disturbance.  Hence  a  disturbance  rejection 
specification  is  determined  based  on  modeled  disturbance  inputs  and  the  basic  QFT  design 
pretense  of  unit  impulse  inputs. 

Since  the  most  severe  disturbance  is  due  to  wind,  see  Appendix  D,  the  disturbance 
specification  is  "tuned"  to  the  wind  input  of  10  ft/sec.  Section  3.4  indicates  a  maximum 
deviation  from  the  nominal  set  position  of  2  feet  in  any  direction  will  confine  the 
receiving  aircraft  to  a  volume  that  will  permit  continued  fuel  delivery.  Therefore  the 
following  disturbance  specification  is  derived.  Given  an  impulse  input  of  magnitude  10 
feet/sec,  the  system  response  will  deviate  no  more  than  2  ft.  Additionally,  the  system 
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will  attenuate  to  half  the  maximum  deviation  in  less  than  1  second.  Equation  (5.1) 
identifies  the  transfer  function  for  the  disturbance  rejection  specification,  and  Figure  5.1 
shows  the  disturbance  rejection  model  response  to  an  10  ft/sec  impulse  input. 

Disturbance  rejection  model  »  - -  (5.1) 

(s  ■*- 1 )  (s  5)3 


Figure  5.1  Disturbance  Rejection  Model  Response  to  10  ft/sec  Impulse 


In  Figure  4.1  the  Lm  of  the  disturbance  rejection  specification  is  superimposed  over 
the  Lm  P(s)  M1SO  loop  plots.  Notice  that  MISO  loops  (2,1),  (3,1),  (3,2),  (1.3),  and 
(2,3)  are  below  the  disturbance  specification  before  compensation  is  applied. 
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5.3  Loop  Shaping 

The  order  of  loop  shaping  is  determined  by  the  amount  of  cross  coupling  each  MISO 
loop  exerts  on  each  other.  As  discussed  in  Section  4.2,  channel  2  couples  strongly  into 
channel  1.  Therefore,  channel  2  is  designed  first,  the  improved  method  is  applied  to 
utilize  the  known  g2  to  recalculate  the  disturbance  bounds  for  channel  1.  Channel  1  is 
then  designed.  Channel  3  is  designed  last  since  it  is  completely  decouple  and  considered 
a  lxl  SISO  system. 

The  bandpass  of  the  plants  are  relatively  low,  a  benefit  of  the  autopilot,  in  shaping 
the  loops  the  overall  system  bandpass  is  designed  to  remain  approximately  equal  to  the 
plant  bandpass.  This  requirement  may  require  tradeoffs  on  meeting  certain  higher 
frequency  bounds. 

5.3.1  Channel  2  Loop  Design 

For  channel  2  plant  case  2  is  chosen  to  be  the  nominal  loop.  Plant  2  is  chosen 
because  through  initial  design  attempts  it  proved  to  be  the  most  difficult  to  shape  around 
the  stability  contour.  A  successful  shaping  of  plant  case  2  guarantees  stability  for  all 
plant  cases.  Templates,  stability  and  disturbance  boundaiies  are  calculated,  see  Appendix 
E.  Composite  bounds  are  formed  in  the  MIMO  QFT  CAD  Package.  The  channel  2 
plants  arc  360  degrees  out  of  phase  between  the  plants  derived  from  aircraft  with  C,  = 
0.2  and  CL  =  0.6.  This  is  evident  from  the  360  degree  wide  templates  and  the  stretching 
of  the  bounds  over  360  degrees.  The  phase  difference  did  not  present  a  problem  as  the 
MIMO  CAD  Package  was  able  to  accommodate  this  scenario.  Compensator  g2  poles  and 


zeros  are  added  to  shape  the  loop.  Channel  2  is  relatively  easy  to  shape  and  proved  to 
have  the  lowest  order  compensator,  g2.  shown  in  Appendix  F. 

Shown  in  Figure  5.2  the  channel  2  loop  easily  satisfies  all  QFT  loop  shaping 
requirements  for  composite  bound  and  stability  contours,  guaranteeing  a  stable  design 
satisfying  the  disturbance  rejection  specification.  Figure  5.3  shows  the  Nichols  plot  for 
all  16  plants.  From  this  figure  the  360  degree  phase  difference  in  the  plant  is  evident. 
Though  there  is  a  phase  difference,  each  plant  correctly  goes  around  the  stability  contour 
indicating  a  stable  design  for  all  plant  cases.  The  following  compensator  is  designed  for 
channel  2. 


(s  +  0.25 )(s  t  0.75 )( s  +  1.2)(s  +  1.3) 
s(s  +  0.98±  l)(s+ 10)(s  +  20)(s  +  120) 
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5.3.2  Channel  1  Loop  Design 

After  g7  is  designed  the  improved  method  is  applied  using  the  equations  derived  in 
Section  2.4.  Utilizing  the  known  structure  of  g2  a  more  accurate  calculation  of  the  cross- 
coupled  disturbance  from  the  compensated  channel  2  to  the  uncompensated  channel  1  is 
achieved.  New  disturbance  and  hence  composite  bounds  are  generated,  shown  in 
Appendix  E.  The  new  bounds  have  smaller  magnitude,  thus,  overdesign  is  reduced. 

For  the  same  reason  as  channel  2,  the  nominal  loop  for  channel  1  is  plant  case  2. 
Again,  as  in  channel  2  the  templates  show  a  360  degree  phase  difference  between  the  two 
plant  cases  of  CL  =  0.2  and  0.6.  But  unlike  channel  2  there  is  a  magnitude  uncertainty 
evident  in  the  channel  1  templates,  see  Appendix  E.  The  magnitude  uncertainty  arises  due 
to  the  strong  coupling  from  channel  2  into  channel  1,  and  also  from  the  difference  in  the 
effect  of  wind  disturbance  between  the  2  classes  of  aircraft  plants  based  on  CL.  The 
plants  of  CL  =  0.6  have  a  larger  wind  induced  disturbance  as  shown  in  Appendix  D. 
Therefore  these  plants  have  not  only  more  external  disturbance,  but  also  larger  cross 
coupling  disturbance. 

The  loop  shaping  is  more  difficult  for  channel  1 .  The  loop  tends  to  curl  at  certain 
frequencies  as  shown  in  Figure  5.4.  The  curling  causes  a  large  change  in  phase  with  little 
or  no  change  in  magnitude.  This  type  of  behavior  makes  it  difficult  to  shape  a  loop  that 
is  stable,  satisfies  the  composite  bound  criteria,  and  maintains  a  low  system  bandpass. 
The  loop  for  channel  1  is  shaped  with  a  compromise  on  the  bandpass.  A  lag-lead 
compensator  is  used  to  "stretch"  the  low  frequency  curl.  Additional  lag-lead  compensators 
are  tried  to  further  "stretch"  the  curl  but  caused  the  loop  to  increase  in  magnitude  as  the 
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frequency  increased.  A  loop  shape  is  finally  achieved  that  satisfies  the  lower  frequency 
bounds,  stability,  and  slightly  increases  the  system  bandpass.  The  channel  1  compensator, 
g„  has  a  higher  order  than  the  channel  2  compensator,  see  Appendix  F.  This  is  an 
indicator  of  the  difficulties  in  achieving  a  loop  shape  that  satisfies  design  criteria. 

The  Nichols  plot  of  all  16  plants  in  Figure  5.5  shows  the  uncertainty  in  the  low 
frequency  range  of  the  plants.  Though  there  is  large  phase  and  magnitude  differences 
between  the  plants  the  QFT  method  is  able  to  achieve  a  design  that  satisfies  stability  and 
disturbance  rejection  for  all  plant  cases.  The  following  compensator  is  designed  for 
channel  1. 


(s  +  0.3)(s  +  0.25±  0.433 )(s  +  3)(s  +  9)(s  +  1.14±  3.747)(s  +  200) 
s(s  +  2)(s  +  0.32±  3.1S4)(s  +  90±  4.02D-o)(s  +  135±  65.38 )(s+  1100) 


(5.3) 
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5.3.3  Channel  3  Loop  Design 


Channel  3  exhibits  none  of  the  channel  1  or  channel  2  characteristics.  The  channel 
3  templates  have  relatively  small  phase  and  magnitude  uncertainty.  There  is  no  coupling 
from  channels  1  or  2  into  channel  3.  The  external  disturbances  have  similar  effects  on 
channel  3  for  all  plant  cases. 

The  lack  of  cross  coupling  disturbance  and  relatively  certain  external  disturbance  is 
evident  in  Figure  5.6  where  the  bounds  collapse  around  the  stability  contour.  The  channel 
3  loop  has  a  tendency  to  curl  up  as  the  frequency  increases.  The  main  difficulty  is  to  add 
compensation  to  shape  the  loop  around  the  bounds  and  stability  contour  at  +180  degrees 
and  then  add  further  compensation  to  keep  the  loop  from  penetrating  the  stability  region 
at  -180  degrees.  To  achieve  stability  the  very  low  frequency  bounds  are  penetrated.  This 
tradeoff  is  considered  acceptable  since  channel  3,  y  position  has  the  largest  margin  of 
disturbance  allowed,  7.5  feet,  as  detailed  in  Section  3.4. 

The  Nichols  plot  of  Figure  5.7  shows  a  very  tight  grouping  of  all  plant  cases.  Again, 
further  evidence  of  relatively  small  uncertainty  in  channel  3.  Notice  the  large  change  in 
phase  with  no  decrease  in  magnitude.  This  is  deemed  acceptable  since  it  occurs  below 
the  zero  dB  line  at  frequencies  below  the  cutoff.  The  following  compensator  is  designed 
for  channel  3. 

(s  +  0.05 )(s  +  0.1  )(s  +•  0.2 )(s  +  0.6)(s  +  1.5±  2.6)(s  +  5)(s  +  30) 

A’3~  s(s  +  2.5D-4)(s  +  0.6±  1.91  )(s  +  10)(s  +  35±  35.707 )(s  +  37.5 ±  64.95) 
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90. 


s-Domain  Open  Loop  Transmissions  for  Channel  3 
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5.4  Closed  Loop  Lm  Plols 

The  overall  equivalent  MISO  system  closed  loop  Lm  plots  are  shown  in  Figure  5.8. 
From  these  plots  you  can  easily  see  that  disturbance  rejection  specification  is  met  for  all 
MISO  loops  except  fo;  as  noted  in  the  low  frequency  portion  of  the  MISO  loop  {3,3}. 


MISO  loop  1 1.1) 
Oj 

-20.  | 


MISO  loop  (2.1) 
O' 

•20.1 


-60 


jSi _ f_(  rad/sec) 

*0.1  1.  10.  100 


Figure  5  8  MISO  Equivalent  System  Lm  Plots 
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The  closed  loop  MISO  plots  of  Figure  5.8  are  an  excellent  indicator  of  success  in  meeting 
the  design  specification. 

5.5  Summary 

In  this  chapter  the  AFCS  is  designed.  Each  loop  shaping  is  detailed,  covering  the 
particular  difficulties  in  shaping  the  loops  for  each  channel.  Also  the  inherent  nature  of 
QFTs  ability  in  handling  large  plant  uncertainties  is  discussed.  Finally  the  Lm  of  closed 
loop  MISO  system  is  shown,  indicating  a  successful  design  given  the  tradeoffs  made.  The 
next  chapter  discusses  the  linear  and  nonlinear  simulations. 
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VI.  Air-to-Air  Refueling  Simulations 


6.)  Introduction 

In  this  chapter  the  compensators  designed  in  the  previous  chapter  are  installed  in  the 
AFCS  and  simulations  are  inn  to  analyze  their  performance.  Linear  simulations  are  run 
for  all  plant  cases  in  MATRIX*.  Nonlinear  simulations  are  for  two  plant  cases,  one  for 
each  CL  =  0.2  and  0.6,  are  performed  in  EASY5x. 

6.2  Linear  Simulations 

Linear  simulation  are  performed  in  MATRIX*  with  the  modeled  external  disturbances 
forcing  the  system  to  deflect  from  the  set  point.  The  simulations  are  executed  in  the 
presence  of  all  external  disturbances  simultaneously. 

The  results  of  the  linear  separation  for  channel  1  (Z  separation)  demonstrate  excellent 
results  with  very  little  perturbation  from  the  set  point.  Figure  6.1  presents  the  channel  1 
response.  The  plots  demonstrate  two  distinct  responses  corresponding  to  aircraft  CL.  The 
aircraft  with  CL  =  0.2  show  a  maximum  perturbation  of  approximately  0.0025  feet.  Also 
the  response  dampens  faster  for  the  aircraft  modeled  with  CL  =  0.2.  The  aircraft  with  C,_ 
=  0  6  deflected  to  a  maximum  value  of  approximately  0.008  feet  with  slower  dampening. 

The  channel  2  (X  separation)  linear  simulation  demonstrates  similar  characteristics 
for  response  based  on  CL.  Again,  excellent  rejection  of  external  dh  turbance  is  achieved 
as  shown  in  Figure  6.2.  CL  =  0.2  aircraft  have  a  maximum  deflection  of  approximately 
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0.025  feet,  while  CL  =  0.6  aircraft  deflect  approximately  0.425  teet  from  the  set  point. 
Recall  that  the  aircraft  with  CL  =  0.6  have  a  larger  uncompensated  perturbation  due  to 
external  wind  disturbance, 

Channel  3  (Y  separation)  has  the  largest  pirturbation  from  the  set  point  in  the  linear 
simulauon,  see  Figure  6.3.  The  maximum  perturbation  in  channel  3  is  approximately  1.9 
feet.  Though  considerably  larger  than  channels  1  and  2,  the  channel  3  perturbation 
remains  within  the  design  specification. 

6.3  Nonlinear  Simulations 

The  nonlinear  simulation  are  performed  in  EASY5x.  EASY5x  has  a  Dryden  wind 
gust  model  preprogrammed  in  the  CAD  package.  The  Dryden  wind  gust  model  is  used 
in  the  nonlinear  simulations  verses  the  disturbance  model  developed  in  Chapter  HI.  Two 
nonlinear  simulations  are  run.  One  representing  an  aircraft  with  CL  =  0.2  and  the  second 
for  CL  =  0.6.  The  nonlinear  simulations  require  considerable  time  to  setup  and  perform, 
therefore,  time  limitation  prevented  performing  a  nonlinear  simulating  for  each  plant  case. 

The  nonlinear  simulations  demonstrate  the  same  excellent  results  that  are  achieved 
in  the  linear  simulation.  The  nonlinear  results  are  consistent  with  the  linear  results,  very 
small  perturbations  for  channels  1  and  2  with  a  larger  deflection  in  channel  3,  shown  in 
Figure  6.4  and  Figure  6.6.  As  in  the  linear  simulations,  the  nonlinear  simulations  are 
within  the  design  specifications. 

Also  presented  in  the  nonlinear  simulation  plots,  Figure  6.5  and  Figure  6,7,  are  the 
control  surface  and  thrust  response  of  the  autopilot.  The  aileron,  rudder,  and  elevator 
responses  are  well  within  the  physical  capability  of  these  devices.  On  the  other  hand  the 
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thrust  requirements  are  probably  beyond  engine  response  capability.  The  engine  response 
is  most  Likely  due  to  the  autopilot  design.  The  autopilot  is  a  "text  book”  design  and  is  not 
very  sophisticated.  A  QFT  design  using  the  actual  C-135B  autopilot  can  probably  achieve 
similar  results  without  extreme  engine  response  requirements. 

0.4  Summary 

The  compensators  designed  in  Chapter  V  are  integrated  into  the  air-to-air  refueling 
AFCS.  Linear  and  nonlinear  simulations  are  performed  with  excellent  results.  The 
system  response  is  within  design  specification.  The  QFT  design  process  worked 
extremely  well  in  designing  the  AFCS  in  the  presence  of  external  output  disturbance.  The 
next  chapter  contains  the  conclusion  and  recommendations  for  this  thesis. 
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Deflection  in  feel 


Figure  6.1  Linear  Simulation  -  Z  Separation  Deflection  all  Plant  Cases 
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X  Deflection  in  feet 
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Figure  6.2  Linear  Simulation  -  X  Position  Deflection  all  Plant  Cases 
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Deflection  in  feet 


Delta  Y  (ft)  Delta  X  (ft)  Delta  Z  (ft) 


Delta  Y  (ft)  Delta  X  (ft)  Delta  Z  (ft) 
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Figure  6.6  Nonlinear  Simulation  -  X,  Y,  Z  Position  Deflection,  Plant  2  CL  =  0  6 
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Figure  6.7  Nonlinear  Simulation  -  Control  Surface  and  Throttle  Response,  Plant 
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VII.  Conclusions  and  Recommendation 


7.1  Discussion 

The  purpose  of  this  thesis  is  to  develop  an  AFCS  to  regulate  the  position  of  a 
cargo/transport  aircraft  relative  to  the  tanker  during  air-to-air  refueling  using  QFT.  To  do 
this  aircraft  models  are  developed  that  envelop  the  structured  uncertainty  of  the  fuel 
receiving  aircraft.  The  uncertainties  are  based  on  dynamic  charges  in  the  aircraft  due  to 
taking  on  fuel  and  the  varying  CL.  Once  the  aircraft  models  are  developed,  an  autopilot 
is  designed  to  control  altitude,  speed,  and  heading  The  autopilot  ensures  a  stable 
minimum-phase  plant  with  a  low  frequency  cutoff.  Disturbance  models  are  developed  to 
identify  the  effects  of  wind  and  refueling  disturbances  on  the  plant.  The  disturbances 
enter  the  system  at  the  plant  output,  therefore,  the  QFT  equations  are  redeveloped  to 
include  the  effect  of  external  output  disturbance.  Since  the  MISO  equivalent  is  not 
diagonal  dominant,  the  improved  method  QFT  equations  are  rederived  to  include  the 
external  disturbance.  The  MIMO  QFT  CAD  package,  developed  by  Mr  Richard  Sating 
is  used  to  design  the  compensators  for  the  control  system.  After  the  loops  are  shaped  the 
QFT  designed  compensator  are  included  in  the  control  system  where  linear  and  nonlinear 
simulations  are  performed.  The  linear  simulation  uses  the  disturbance  model  to  generate 
the  external  output  disturbance.  The  Dryden  wind  gust  model  is  used  in  the  nonlinear 
simulations. 
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7.2  Conclusions 


M1MO  QFT  CAD  Package  An  outstanding  QFT  design  tool.  The  CAD  package 
greatly  enhanced  the  design  process.  Though  problems  were  initially  encountered, 
they  were  quickly  solved  and  the  package  proved  to  be  very  robust  and  easy  to  use. 
QFT  External  Output  Disturbance  Equations  The  successful  design  of  compensators 
based  on  these  equations  and  the  subsequent  simulations  validate  that  the  equations 
are  correct.  These  equation  are  now  part  of  the  MEMO  QFT  CAD  Package  and 
should  prove  to  be  valuable  in  future  QFT  design  work. 

Air-to-Air  Refueling  AFCS  The  excellent  results  from  the  linear  and  nonlinear 
simulations  demonstrate  the  very  real  possibility  of  being  able  to  implement  this 
system  on  Air  Force  aircraft. 

7.3  Recommendations 

The  successful  simulations  the  QFT  designed  air-to-air  refueling  prove  the  technology 
exists  to  develop  this  AFCS.  Further  research  should  be  accomplished  using  more 
modern  aircraft  with  actual  inservice  autopilot.  Upon  successful  conclusion  of  this 
research  ,  the  Air  Force  should  give  serious  consideration  to  implementing  this  system. 
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Appendix  A.  C-135  Nondimensional  Stability  Derivatives 

This  appendix  contains  the  nondimesioanl  stability  derivatives  used  to  generate  the 
C-135B  bate  aircraft  model. 

A.l  Nondimensional  Stability  Derivative  Definitions 

A.  1.1  Lonitudinal 

C„  -  Aircraft  drag  co^. .  ..;ient 
C,  -  Aircraft  lift  coefficient 
C,„  -  Lift  curve  slope 

C(A.  *  "Change"  in  lift  due  to  a  change  in  eievator  deflection 

Cm;:  -  "Change"  in  pitching  moment  due  to  a  change  in  alpha 

Cm„  *  "Change"  in  pitching  moment  due  to  a  change  in  alpha  rate 

Cmq  -  "Change"  in  pitching  moment  due  to  a  change  in  pitch  rate 

CmfH.  -  "Change"  in  pitching  moment  due  to  a  change  in  elevator  deflection 

A. 1.2  Lonitudinal 

C>(!  -  "Change"  in  side  force  due  to  a  change  in  beta 

C>p  -  "Change"  in  side  force  due  to  a  change  in  roll  rate 

Cyr  -  "Change"  in  side  force  due  to  a  change  in  yaw  rate 

C)flI  ‘  Change"  in  side  force  due  to  a  change  in  rudder  deflection 

C>fl,  -  "Change"  in  side  force  due  to  a  change  in  aileron  deflection 

Cl(,  -  "Change"  in  rolling  moment  due  to  a  change  in  beta 
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C|P  -  "Charge"  in  rolling  moment  due  to  a  change  in  roll  rate 

Cb  -  "Change"  in  rolling  moment  due  to  a  change  in  yaw  rate 

Cl6,  -  "Change"  in  rolling  moment  due  to  a  change  in  rudder  deflection 

CI6J  -  "Change"  in  rolling  moment  due  to  a  change  in  aileron  deflection 

CnP  *  "Change"  in  yawing  moment  due  to  a  change  in  beta 

Cnp  -  "Change"  in  yawing  moment  due  to  a  change  in  roll  rate 

Q,  -  "Change"  in  yawing  moment  due  to  a  change  in  yaw  rate 

Co6r  -  "Change"  in  yawing  moment  due  to  a  change  in  rudder  deflection 

Cn6ll  -  "Change"  in  yawing  moment  due  to  a  change  in  aileron  deflection 


Stability  Derivative 

Value 

c>P 

-0.74485 

Cvr 

0.410 

0.253 

C>6  a 

-C.0074 

C* 

-0.3550 

C|6, 

0.0160 

C„6r 

-0.109 

C„6, 

0.00160 

-4.641 

C,„ 

-0.1977 

-0.8309 

-0.67609 

Table  A.l  Constant  Stability  Derivatives  for  all  Plant  Cases 
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Stability  Derivative 


Value 


CL 

0.2 

0.6 

-0.237 

-0.196 

-0.17762 

-0.22059 

c„ 

0.108 

0.192 

Qftr 

0.03724 

0.0258 

Qifl 

0.15183 

0.16043 

Cnp 

0.0 

-0.006 

Cw 

-0.183 

-0.187 

Q) 

-0.014 

-0.03 

Table  A.2  Stability  Derivatives  that  Change  with  CL 


HI 

wmwMWm, 

Center  of 
Gravity 
(%  MAC) 

Qnd 

c 

(xlO6) 

(xlO6) 

h 

(xlO6) 

1 

160.166 

28.8 

-6  10 

-M.25 

2.0 

2.70 

4.575 

2 

207.316 

28.0 

-6.13 

-14.5 

2.725 

2.70 

5.375 

3 

210.189 

30.1 

-6.03 

-14.275 

3.15 

2.75 

5.85 

D 

245.500 

22.9 

-6.35 

-15.025 

3.175 

2.855 

6.0 

5 

253.500 

26.0 

-6.22 

-14.725 

3.1875 

3.0 

6.1 

6 

263.500 

23.8 

-6.30 

-’4.925 

3.2 

3.05 

6.19 

7 

275.500 

24.1 

-6.30 

-14.915 

3.525 

3.10 

6.5 

8 

277.500 

24  8 

-6.27 

-14.835 

3,55 

3.105 

6,51 

Table  A. 3  Stability  Derivatives  that  Change  with  Weight  and  Center  of  Gravity 


A-3 


Appendix  B.  Plant  Transfer  Functions 


This  appendix  presents  transfer  functions  for  all  plants  developed  with  EASY5  and 
MATRIX,. 

B.l  Plant  Case  1  -  CL  =  0.2  Gross  Weight  =  160,666  pounds 

a  2.05P-3Q*  0,286 1 4.346)0  *  17023)0*  43.763  ±  73.228)0 *  159.108)0*  499.10*)  _ 

is-  2.068  ±6.221)0 *  0.067)0*  0.754 *  1.208)0*  52394)0*  10169)0*  70040±7t.419)0*  119.979) 

b  2.42D  *  1 3Q  *  1 7.223)0  *  43.763 ± 73.228)0  *  159. 108)0  ♦  249.323 ± 433.897)0  -  499. 1 04) 

Qn  '  sis -621 D  *4  i  3.640  *  4)0  *0.001)0  *  1)0  *  50)0  *  444  .55 1 7080  *  4)0  ♦  6.350  *4 1 3.640 .4) 

<7u*  0 

b  _ -12.732Q*  0.269  ±4,346) 

4ll‘  0-  0.023  ±0.013)0*  0.977  ±1.7 18)0*  10.031) 

.  49.866Q  *  0.289  ±4.346) _ 

Qa  iO*  0.099)0*  0289  ±4.336)0*  501^11) 

9«-  0 


9±i*  0 


1nm  0 


.  2.296  0  -  0  037  1 1 .739)0  *  0-2)0  *  0-289)0  ♦  5.620  ±  6. 1 10) _ 

*  sis  * 0  1 46  ±  0072)0  -  02185)0  *1  1 15)(0  * 5.430  ±  5.95 1  )0  *  7. 1 36)) 


B. 2  Plant  Case  2  -  CL  =  0.6  Gross  Weight  =  160,666  pounds 


2.780-30  *  0471  ±4.346)0*  17.223)0*  43.123  ±  72.056)0  *  W  600)0*  369.868) 
q"'  if  0.268  i  19.368)0  *  0.C22  ±0.820)0  *  0.067,0  •  1. 362)0 ”  14.244)0*  70.799  ±  71  .556)0-  119674) 
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2. 1 830  - 1 6 (f  -  0.47 1  ±  4.364)(/  *  1 7298)0  *  43. 1 23  ±  72056)C»  *  l  \  6*-399)ft  *  249-279  *  433.59  i  )(j  *  369.868) 
iC»  -  0.00*)(/  *  0286 *  4.340)(j  *  1  )C*  ♦  50) 


«  u* 


9 o’  0 


.  _ -11853(4*  0.477  *4.36*) _ 

Ci  0.025  *0.040)(s 4  0.759*  1.818)C»*  10.029) 


.  37.163(4*  0.477*4.364) _ 

4(j*  O.IOl)Cs *  a473*4.J41)(4*  376574) 


In"  0 


9ji*  0 


0 


,  _ 22 15  (j- 0.095 * 1 .763  )(4*02)(4*0.498)(4  *5.615  *  6391) _ 

4 CS  *  0. 1 38)0  -  a  1 56)(i  *  0.344)0  - 1  .444)0  *  1 .957)(i  *  5.392  *  5.795 )(*  ♦  7.01 9) 


5.3  Plant  Case  3  -  CL  =  0,2  Gross  Weight  =  207,316  pounds 


.  2.04D  -  3(3*  0-208*  3.697X4 ♦  172275)0 *  43.356*  72.337)0*  Ki2.161)U*  413.414) 

Qu  (i-  2.119*6292)(j  *  0.062)(j ♦  0.704  *  0891  )0»*  5.40)0*  10.l78)(i*  70.032  *  71.430)0*  119.994) 

.  1-3080  - 160  *  17.275)Cs  *43.356  *  72537X4  ♦  162. 161)(t  *249  277  *  433.702X?  *  413.414) 

,u‘  jC»-6.Q02)(j*1)Cj*50) 


9o*  0 


_  _ -13.62/ fj*  0209*3.698) _ 

<fj|"  (*  0.029)(i-  0.019) i*  ♦  0.898*1.701)0*  10030) 


_  _ 41.405(j  *  0.209*3.698)  _ 

j(j*  0.099)C»  *  0.21 1  *  3090)(4  •  418.121) 


9u*  0 
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9)1  m  0 


9nm  0 

.  2.1 13  Q- 0.060  ±  1.483)0  *02)0*5.445  *4.568) 

9u  ~  Mi j  ♦  0.154 1  0.069)0  ♦  1275  ±  0.444XJ  *5210 ± 4.318)0  ♦  7.739) 

BA  Plant  Case  4  -  CL  -  0.6  Gross  Weight  =  207,316  pounds 

m  _ 2.6QD-3Q*  P.380*  3.718)0*  17,330)0*  42.804*71.411)0*  168.723)0*  322.491) 

9"‘  is-  0.748  ±  19295)0  *  0.062)0 ♦  0.076 1 0.684)0  *  1.263)0*  14283)0*  70.628 1.71.744)0*  119.960' 

.  -349.466Q  *  0.381 1 3.7 1 8)0  ♦  1 7.33 1  )Q  *  42.806  * 7 1 ,4 1 1  )p  *  1 68.724)0  *  249.302  *  433.469)0  *  322.4906) 
’  /fit-  1.695D  *9)0-0.007)0  *0207  ±  3.690)0  *  1X»  *50)0  *  1-69523  *9) 

9u*  0 

.  _ -12.7380  *  0.384 ±3.7 18) _ 

9il'  0-  0.027  ±0.032)0  *  0.611  ±1.815)0*  10.029) 

.  _ 32.740Q*  0,384  ±3,718) 

JO*  0.100)0*  0.383 ±  3.702)0  *  332.784) 

9a*  0 

9),m  0 

9jj*  0 

.  _ 2Q52Q -0.133  ±  1.530)0  *  0.2)0  *  0.5 18)0  *  5.422  ±  4,802) 

"  JO  *  0. 126)0  *  0 176)0  *  0.347)0  *  1248  ±  0625)0  *  5228  ±  4.108)0  *  7.636) 
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B.5  Plant  Case  5  -  CL  =  0.2  Gross  Weight  =  210,189  pounds 


_ 2,010-30+  1207*3674)0  +  17480 )Q  +  43.321*72.474)0  +  162,466)0 »  407.121) _ 

O-  11223  *  6052)0+  0.062)0+  0.702*  0.875)0  ♦  5.368)0+  10.172)0+  71031*71.429)0+  119.994) 


_  8,8540+40  +  17080)0  +  43421  *71,474)0  +  162.466)0  +  249476*433686)0  +  407.121) 
,u  "  *0  -4.1850  +6*74490  +6)0+0:002)0  ♦  1)0  +  50)0  ♦  8.3710  +9) 


?o"  0 


_ -13.365Q+  0409*3.674) _ 

O'  0.030)0+  1018)0  +  0.874*1689)0  +  10.030) 


.  40.792Q  +  0409*3674) 

40+  0.099)0+  1210*3467)0  +  412.067) 


«oj-  0 


«3i*  0 


«s*  0 


.  1 .833  Q  -  0.079  *  1 6 12)Q  + 04)0  +  5-406  *  3,822) 

40  ‘  0.155*1069)0  +  1-143  *  0571)0  +5410*3.523)0  +  7.981) 


B .6  Plant  Case  6  •  CL  =  0.6  Gross  Weight  =  210,189  pounds 

_  1298Q+  0.384*3695) _ 

O*  0.062)0+  0.080*0.679)0+  1460)0+  U  4.471) 

.  61750  -  30  +  0.384  *  3695) 

,u"  40  '0.007)0  +0407 *3.667)0  *  1) 

«ji*  0 


_ -13.475Q+  0484*3694) _ 

O-  0.027+0.032)0+  0489*1.799)0+  10.028) 
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s  _ 32^290  +  0.384  ± 3.695) _ 

4a  jO  *  0. 100)0  +  0.382  ±  3.678)0  +  327.733) 


?»-  0 


9ji“  0 


9r*  0 


_ 1.776Q  -0.150  i  1.338)0  ♦  02)0  +0324)0  +-5.374  i  4,070) 

jO  *  0.1 25)0  *  a  1 77)0  *  0.349)0  *  1 2 13  ± 0.720)0  *  3. 1 69  ±  3^62)0  *  7.878) 


B.7  Plant  Case  7  •  CL  =  0.2  Gross  Weight  =  245,500  pounds 


.  _  3.232Q*  0.170.3,361) _ 

4,1  ’  o*  0.058)0+  0.677 ±0.710*0+  8.867)0+  119.936) 


.  3.244X)-3p +■0.170  ±3.361) 

4u*  JO -0.003)0‘0.169±  3.351)0  +  1) 

9ji*  0 


.  _  -13.193Q+  0.170±3.361) _ 

4il  ’  O-  0.037)0-  0.013)0+  0.840 ±1.639)0+  10.029) 


9a* 


_ 35.848Q  +  0.170  ±3.361) _ 

iO+  0.098)0+  a  172  ±  3.334)0  +  363.322) 


9u*  0 

9ji  "  0 


9r*  0 


_ 1.983Q  -0.07 1  ±  1 ,346)0  +  02)0  +0266)0  +  3.381  ±  3.799) 

JO +0.138)0 ‘0.067)0 ‘0^73)0  + 1.099  ±0591)0  ‘  3^07  *  3,486)0  +  8.049) 
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B8  Plant  Case  8  *  CL  =  0.6,  Gross  Weight  =  245,500  pounds 


. _ 0.294Q*-.-  0.341 1 3.382) _ 

11  <**  0.058)0**  0.106  ±0.606)0**  1.213)0**  128.756) 


.  5.285D- 3  (**0.341  ±3.382) 

,u’  *0* -0.008)0 *O.169±  3.352)0*  *  1) 


flu”  0 


_ - 1 3.296 Qt*  0.341  ±3.382) _ 

(*-  0.027 ± 0.027)0* *  0.5 17  ±  1.575)0**  10.027) 


fla* 


28.999 (**  0.341  ±3.382) 
s(s  *  0.100)0**  0.339  ±3.367)0*  ♦  295.845) 


flu*  0 


flji*  0 


flu*  0 


_ 1.924Q* -0,141  ±  1.4Q7)C**0.2)0**QS30)(**5.352  ±4.014) 

M  *(*  *0.123)0*  * 0 184)(*  *0.349)0*  *  1.163  ±0.735)0*  *5.165  ±  3i33)(*  *  7.957) 


5.9  Plant  Case  9  -  CL  =  0.2,  Gross  Weight  =  253,500  pounds 


.  _ 4.984  Q*  *  0,168  ±3.313) _ 

9“  C*  *  0.058)0*  *  0.672  1  0.679)0*  *  8.89S)(*  *  1 1 7.665) 

«  5.0780-30**0.168  ±3.313) 

9u  *(*-0.003)0*  *0.167  ±3.303)0**  1) 


flu*  0 


.  _ -I2.S19Q**  0.168 ±3.313/ _ 

9n  0*  0.038)0*-  0.012)0**  0.793  ±1603)0**  10.027) 
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In  * 


_ 34,607(5*  0.168  ±3.313) _ 

5(5*  0.098)(s*  0.170±3.306)(5*  331.087) 


9u*  0 


In'  0 


9n*  0 


1  .973(5  -  0.070 1 1 .32 1  )(5  *02)(5  *  026t)(5  *  3.379 1 3.780) 

5(5  *  0. 1 59  ±  Ol066)Cj  *  0  272)(y  *  1,073  ± 0.606)(5  *  3 .203 1 3.467)(j  *  8.094) 


B.10  Plant  Case  10  -  CL  =  0.6,  Gross  Weight  =  253,500  pounds 


_  0.279(5*  0.339 1  3.334) _ _ 

q"‘  (s*  0.038)(5*  0.111  iO-59l)U*  1203)(j ♦  126J02) 

.  3. 1180 -3(5  *  0.339  ±3.334) 

5(5 -0.008)U *0.167  ±3.304)(5*1) 


9  13*  0 


_ -12617(5*  0.339*3.334) _ 

(5  -  0.028  ±  0.026)(5  *  0.476  ±  1.713)(5  *  10.026) 


,  27,906(5*  0.339  ±  3.334) 

5(5*  0.100)(5*  0.337 ±  3.319)(5  *  285.044) 

9u“  0 


95.*  0 


95J"  0 


%  _ 1.913(5-0.133»1.3t'3)(j*a2)(5*0J31)(5*5.349±3.992) _ 

*  5(5  *0. 1 22)(5  *  a  183 )(5  *  0.349)C5  *  1.1 35  1  0.74 6)tf  *  5. 164  ±  3-220)0  *  8.007) 
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B.12  Plant  Case  12  -  CL  =  0.6,  Gross  Weight  =  263,500  pounds 

_  _ 02760*  0.330  ±3264) _ 

(i*  0.037)0*  0.117  ±0374)0*  1.193)0*  129.283) 

.  4.9240-30  *0.330  ±3.264) 

,u  JO  -0.008)0  *0.1 38  ±  3.234)0  *  1) 

?u*  0 

.  _ -12438Q*  0.330  ±3264) _ 

,J1  O- 0.028  ±0.023)0*  0.438  ±  1696)0*  10.023) 
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_  _ 27.074(f«  0,330*3264) 

<J“*  j(j*  ai00)(i*  0.325 1 3250)(s *  276.840) 


flu*  0 


flu*  0 


flu*  0 


_ 1.916(5  -0.135 1 1.356)(»  *  0.2)0  *0.532)(5  *  5.346  *  3.971) 

5(5 * 0. 1 22)(5  *  0 1 87)(i * 0.349)(j *  1.1 11  ± 0.7 J4)(5 *5.162*  3206)0 * 8.050) 


B.13  Plant  Case  13  -  CL  =  0,2,  Gross  Weight  =  275,500  pounds 


_ 4,870(5  *  0,153  ±3.164) _ 

(5*  0  056)0  *  0.659  ±  0.60 1)(»*  8.909)0*  124.181) 


_  4.6720  -  3(5  *  0.163*3.164) 

Vu  '  5(5  0.003 )(.*  *0.152 1 3.153)(5  *  1) 


flu*  0 


. _ -12.232C**  0.153*3.164) _ 

’  C5-  0.040)(5-  0.010)0*  0.763. 1569)(j*  10.027) 


fl«‘ 


_ 32.169(5*  0.153 1 3243) _ 

5(5  *  0.098)(5*  ai55.3.157)(5  *  327.051) 


flu*  0 


fl). 


0 


0 


1.875Q- 0.072. 1266)0*02)0  *0258)0  *5.347  ±3248) 

*  5(5  *0.162  iO065)(5  *0269)(5  *0.990  ±  0.647)0  *  5.1 69  ±  3.87 1 Xj  *  8.240) 
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B.14  Plant  Case  14  -  CL  =  0.6,  Gross  Weight  =  275,500  pounds 


fn* 


_ _  0.2730*  0.321  ±3.185) _ 

0*  0.056)0*  0.123 ±0J55)O*  1.182)0*  133.002) 


4.7080-30*0.321  ±3.185) 
jO  0.009)0  *  0.151  ±3.154)0*1) 


«u*  0 


.  _ -12.327Q  *  0,321  ±3.185) _ 

*JI"  0-  0.028  ±0.023)0  *  0.429 ±  1.680)0  *  10.025) 

.  _ 26.194Q*  0.321  ±3.185) _ 

i(j«  0.100)0*  6. 320±  3.171)0  *  268.140) 


9um  0 


9i,‘0 


0 


_ 1.820Q  -0.163±  1.333)0  *  0.2)0  *0.539)0  *5.31 3  ±  3.762) 

*  JO  *  0. 1 2 1  )0  *  0. 1 89)0  *  0.35 1 )0  *  1 .048  ±  0.779)0  *  5. 1 28  ±  2J73)0  *  8. 155) 


B.15  Plant  Case  15  -  CL  =  0.2,  Gross  Weight  =  277,500  pounds 

.  4.866Q *  0.152±3.151) _ 

?n  O*  0.056)0*  0.658  ±0-594)0*  8.910)0*  124.910) 


4.6390  -  3Q  *0.152±3.151) 
jO -0.003)0  *  0.151  ±3.140)0*1) 


flu*  0 


_ -12.224Q*  0l}2±3.151) _ 

0-  0041)0-  0.010)0*  0.757 ±1.568)0*  10.026) 
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_  _ 31,972(1*  0.152  1  3.151) _ 

,a"  *(i'  0.098)(i*  0.154i3.144)(j  +  325.094) 

«»■  0 


0iim  0 


On  -  0 


.  1-879& -0.G71  *1.259)(J*0.2)(i*  0.257)(i  *5.344  1 3.210) 

'  1(1  -0.162  ±0.065)(i  * 0-269)C l  * 0.988 t 0.648)(i *  5 . 1 67  ±  Z 826)(i  *  8.242) 


B.16  Plant  Case  16  -  CL  -  0.6,  Gross  Weight  =  277,500  pounds 


0,273(1*  0.321  ±3.172) 

(f*  Q.0S6)(i *  0.124 ±0-552)(i *  1.180)(j*  133.736) 


,  4.6750  -3(i  *0,321  1 3.172) 

Qa'  1(1 -0.009)(s  *0.151  }?.U2)(5*1) 


On-  0 


.  _ -1Z319(i*  0.32U3.172) _ 

'  is  -  0  028  1  Q.023)(s  *  0.421  ±  1.680){j  *  10.025) 


. _ 26.060(1*  0.321  ±3.172) 

J(l*  0.099)(i  •  0.319 1  3.1J8)(i *  266.819) 


0 


0s,‘  0 


0n‘  0 


_ 1.824(i  -0.163 1 1.327)0  *  0.2)(i  *k539)(i  * 5.310  1  3.420) 

j(i  *0.12l)(i  *ai89)(i  *0.351)Ci  *  1.046 1 0.780)(i  *  5.126  ±  Z520)y  *  8.157) 
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Appendix  C.  C-135B  and  Autopilot  Input  Response 


This  appendix  contains  the  time  response  plots  for  the  C-135B  and  autopilot. 
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Figure  C  6  Rudder  Control  Response  to  1  deg  Step  Input  -  CL  =  0.6 
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Figure  C.8  Aileron  Control  Response  to  1  deg  Step  Input  -  CL  =  0.6 


Appendix  D.  C-135B  and  Autopilot  Disturbance  Response 


This  appendix  contains  the  time  response  plots  for  the  C-135B  and  autopilot  in  the 
presence  of  external  disturbance. 
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Figure  D.3  X,  Y,  Z  Response  to  Lateral  Wind  Disturbance  -  CL  =  0.2 
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Figure  D.4  X,  Y,  Z  Response  to  Lateral  Wind  Disturbance  *  CL  =  0.6 
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Appendix  E.  Templates  and  Boundary  Plots 


This  appendix  contains  the  template  plots.  It  also  presents  the  stability  and 
disturbance  boundaries  used  to  form  the  composite  bounds. 
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Templates  -  Channel  2  -  Nominal  Plant  Case  2 
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Appendix  F.  QFT  Compensators 


This  appendix  contains  the  compensators  designed  in  Chapter  V. 


F  I  Channel  1  compensator,  g] 

(s  +  0.3)(s  +  0.25 ±  0.433 )(s  +  3)(s  +  9)(s  +  1.14±  3.747)(s  +  200) 
s(s  +  2)(s  +  0.32 ±  3.184)(s  +  90±  4.02D-6)(s+  135±  65.38)(s  +  1100) 


F.2  Channel  2  compensator,  g2 

(s  +  0.25  )(s  +  0.75)(s+1.2)(s+  1.3) 
s(s  +  0.98 ±  l)(s+ 10)(s  +  20)(s+ 120) 


F.3  Channel  3  compensator 

(s  +  0.05 )(s  +  0.1  )(s  +  0.2 )(s  +  0.6)(s  +  1.5±  2  6)(s  +  5)(s  +  30) 
s (s  +  2.5D-4)(s  +  0.6±  1.91  )(s  +  10)(s  +  35±  35.707 )(s  +  37.5±  64.95) 
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